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We propose a minimal model resolving a puzzle of enigmatic correlations observed in sodium- 
rich NaxCo02 where one expects a simple, free motion of the dilute S — 1/2 holes doped into 
a band insulator NaCo02. The model also predicts singlet superconductivity at experimentally 
observed compositions. The model is based on a key property of cobalt oxides - the spin-state 
quasidegeneracy of CoOe octahedral complex - leading to an unusual physics of, e.g.^ LaCoOs- We 
show that correlated hopping between t2g and Cg states leads to the spin-polaron physics at x ^ 1, 
and to an extended s-wave pairing at larger doping when coherent fermionic bands are formed. 

PACS numbers: 71.27.+a, 74.20.Mn, 74.70.-b 



I. INTRODUCTION 

Recent studies boosted by the discovery of water- 
induced superconductivity (SC) in NaajCo02^ ex- 
posed many remarkable properties of these compounds^ 
such as a spin-sensitive thermopower^, unusual charge 
and spin orderings^i^i^i^, very narrow quasiparticle 
bands^i^iiii^i^ etc. While strongly correlated nature 
of Naa;Co02 is no longer at doubt, the mechanisms 
by which the correlated electrons design such an exotic 
phase diagram^ are not fully understood even on a qual- 
itative level. 

Superconductivity of cobaltates has low Tc — 5 K. 
However, the identification of the pairing mechanism is a 
problem of principal importance, because this may shed 
light on the other puzzles of Naa;Co02 as well. Moreover, 
hopeful comparisons with the high-Tc cuprates have been 
madeiii^, noticing that NaajCo02 consists of C0O2 layers 
with 5* = 1/2 Co'^^ ions doped by 5 = Co^+ charge 
carriers, an apparent t2^-band analog of the cuprates. 
A triangular lattice formed by Co ions, providing favor- 
able conditions for a realization of the resonating- valence- 
bond (RVB) idea^, has been also emphasized. 

However, it was quickly realized that: {%) the phase 
diagram of Naa;Co02^ is radically different from that of 
cuprates; {%%) SC dome is located at valence composi- 
tions closer to Co^+(5' = 0) rather than Co'^+(5' = 1/2) 
(z.e., at average valences ^ 3.50) not favorable for 
RVB-theories^^. Further, Naa:Co02 shows magnetic or- 
der at X > 0.75 (besides a particular one at x = 0.5^^^) 
which is counter-intuitive because density of Co^+ spins 
oc (1 — x) decreases at large x. These observations make 
it clear that the origin and functionality of strong corre- 
lations in cobaltates and cuprates are very different. 

In this paper, we propose a model for strong corre- 
lations that operate over the entire phase diagram of 
Naa;Co02 and lead to SC optimized near the valency 3.4 
as observed. First, we consider a single hole doped in 
NaCo02 and show why its behavior is radically different 
from that of a free carrier embedded in a band insulator. 
Considering then a Fermi-liquid regime of NaajCo02, we 



demonstrate how an unusual, kinetic energy driven pair- 
ing emerges in the model. 

The model is based on the following points (none is 
present in cuprates): {%) typically, Co^+ ions in the 
octahedral environment possess also low-lying magnetic 
states, e.g. t^gC^g = 1 or t^gC^g 5 = 2; (n) in the C0O2 
planes with 90° Co-O-Co bonds, the correlated 5* = 1 
spin states are strongly coupled to the ground state via 
the intersite t2g ^ eg hopping (see Fig. [1])^^. In other 
words, the magnetic configuration of Co^+ ions is acti- 
vated once the mobile Co^+ holes are added in NaCo02. 
A dynamical generation of ^2c/^^ S = 1 states by a hole 
motion converts it into a many-body correlated object - 
the spin-polaron. At larger density of Co^+, we elimi- 
nate a virtual 5* = 1 states perturbatively, and find an 
effective model in a form of spin-selective pair hopping 
of electrons. The correlated hopping energy is optimized 
when holes are paired and condense into a SC state. 

Spin-state quasidegeneracy of cobalt ions is well 
known, LaCoOs being a textbook example^!. Accord- 
ing to Ref. [22|, magnetic states are in the range of 
- 200 - 400 meV {S = 1) and > 50 meV {S = 2) 
above the ^2^ 6* = ground state (without lattice re- 
laxations). A balance between the crystal- field, Hund's 
coupling and pd-covalency is easily tuned and latent mag- 
netism of Co^+ living in virtual states can be activated, 
e.g., by nonmagnetic dopin g^^^^^ . In oxides with 180° 
d-p-d bonding as in LaCoOa, this process leads typically 
to a ferromagnetic metal stabilized by an electron pro- 
moted into broad eg bands^^. New element of Naa;Co02 
is the 90° d-p-d bonding where the eg-eg hopping is sup- 
pressed. Instead, a large overlap between the neighbor- 
ing eg and t2g orbit als is dominant. A curious situation 
which arises is that while Co^+ ions are nonmagnetic in 
NaCo02, their S = 1 ^2^6^ configurations are dynam- 
ically generated in a doped case by the strong t2g-eg 
hopping. 

A resume is that a low- lying magnetic states of Co^+, 
accessible for electrons via the intersite hopping, provide 
an extra dimension in physics of Naa^Co02. In Sec. II, 
we design a model incorporating this idea. Based on this 
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model, we demonstrate in Sec. Ill that a hole doped into 
the band insulator NaCo02 behaves in fact as a magnetic 
polaron dressed by the spin-state fluctuations of Co^+ 
ions that are excited by hole motion. Sec. IV derives the 
interaction between holes, mediated by virtual spin-state 
excitations of Co^+ ions, in a Fermi-liquid regime at fi- 
nite hole densities. We also discuss there the relevance 
of these interactions to the spin ordering, and find signa- 
tures of 2/ci?-instabilities. Finally, we focus in Sec. V on 
the superconductivity and discuss symmetry and doping 
dependencies of pairing instabilities within our model. 
Sec. VI concludes the paper. 




a-bond: t (yz ^ 3x^-r^) 



b-bond: t {zx ^ 3y^-r^) 
c-bond: t {xy ^ 3z^-r^) 



II. MODEL HAMILTONIAN 



The t2g orbitals in NaajCo02 split into aig 



{dxy + dyz + dzx)/V^ and e' 



e^^'-^dzx)/V^ states {(p = 27r/3). The photoemission 
experiments^^ ^^'^^^^^ show that a single band, derived 
mostly from the aig orbitals, is active near the Fermi 
level (see Ref. |25| for the orbital-selection mechanism). 
Therefore, we base our model on the aig = f hole states 
(its three-band version will be presented elsewhere^^). 
Valence fluctuations d^d^ d^df within the low-spin 

t2g manifold read then as Ht = -t^ijcr f]crficr^ where 
t = 2^0 /3 and to = tj^t-j^/Apd is the overlap between 
t2g orbitals^^ (hereafter, a hole representation is used). 
Our crucial observation is that the t2g-eg hopping i = 
tcrt-j^/Apd^ which uses a stronger a-bonding path with 
t(j/t-j^ ~ 2, leads to more effective valence fluctuations. 
The hopping geometry is depicted in Fig. [TJ The near- 
est neighbor (NN) Co ions and two O ions binding them 
determine a plane which could be labeled a, 6 or c ac- 
cording to the Co-Co bond direction. With respect to 
this plane, the t-hopping couples the in-plane t2g orbital 
to the out-of-plane Cg orbital. 

The i process generates S = 1 state of Co^+ com- 
posed of a t2g hole and an Cg electron; we represent 
it by T operator (low-spin = ^2^eJ state is much 
higher in energy and can be ignored^^). T is speci- 
fied by its spin projection and the Cg orbital 7 created 



by i hopping, i.e., Tl^^^ 
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show our minimal model for Naa;Co02: H^_^ 
where Ht is as given above, while 
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FIG. 1: (a) Electron hopping from t2g to Cg orbital via oxygen 
atoms in the case of 90° bonds creating S — 1 t^g^g configu- 
ration of Co^^. The t2g orbital laying in the plane (here the 
xy orbital) couples to the out-of plane Cg (3z^ — r^) orbital, 
(b) The coupling to the planar orbital is zero because of the 
destructive interference of two channels, (c) Bond directions 
and corresponding angles in the hexagonal lattice of Co ions 
and t-active orbitals on these bonds. 



"picks-up" the spin correlations in virtual states. The in- 
dex 7 is determined by the orientation of the {ij) bond 
according to the rules in Fig. [U^c). The overlap between 



orbitals specified by 7 and 7' is (7I7') = cos( 
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Consequently, the excitations inherit the same over- 
lap: {T^T^,) (X (7I7O. The T-excitation energy Et 
is determined by all the many-body interactions within 
the CoOe complex (Hund's coupling, p — d covalency, 
crystal field, . . . )^. This is a free parameter of the 
model. Experimentally, 5* = 1 states of CoOe complex 
in perovskite compound LaCoOa are found at energies 
Et ~ 0.2 — 0.4eV^^ as already mentioned in the introduc- 
tion. Based on this observation, we will use in this paper 
a representative value Et — 0.3 eV for layered cobal- 
tates. In units of aig hopping integral t 0.1 eV (which 
follows from the band structure fit to 0.15 eV— ), this 
translates into Et /t = 3 adopted below in our numerical 
data. (In principle, we expect some material dependence 
of Et as it is decided by the balance of several competing 
interactions. It is therefore highly desirable to quantify 
a multiplet structure of CoOe-complex in Naa^Co02 as 
done in LaCoOs^^). For the ratio of the hopping ampli- 
tudes t and to, we set t/to = 2 as t^/t 
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moves an electron from Co^^ to Co^^ - producing a 
t2g hole on site j - and replaces the t2g hole on site i by 
a complex excitation T. Making use of the t2g-eg hop- 
ping (the largest one for 90° Co-O-Co bonds), an electron 



III. SPIN-STATE POLARON 

It is instructive to consider first a single hole doped in 
NaCo02. With Ht alone, it is just a usual plane- wave 
having nothing common with what is actually seen in 
Naa;Co02 at large x. Things change radically when the 
H^ is switched on: now, a hole generates a multiple of 
T excitations [see Fig.[2fb)], and spin-polaron physics of 
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FIG. 2: (a) Spectral functions of a Co^^ hole doped in 
NaCo02 along M-F-K path in the Brillouin zone [inset (ii)]. 
Bare and renormalized dispersions, measured both from the 
chemical potential for better comparison, are shown. (The 
polaron binding-energy shift is Eb ~ 2.6t). (i) The diagram 
describing the dressing of a hole by 5* = 1 T-excitations. 
(b) Density of states compared to that of the bare band. Vir- 
tual processes associated with the polaronic spectral features 
are sketched: every use of t hopping channel generates S = 1 
states of Co^^ behind the hole. 



a typical Mott insulator emerges. Given that NaCo02 
itself is nonmagnetic at all^^, the correlated behavior of 
doped holes in layered cobaltates has been a mystery; it is 
resolved here by invoking a "virtual Mottness" of cobalt 
oxides2S hidden in their low- lying magnetic states. 

Polaron physics is evident from the spectral functions 
in Fig. [21 We have employed a self-consistent Born ap- 
proximation for the selfenergy, which then takes the form: 



k 



Er-^k- - Et) + 



(2) 



Here, F^ = ^ + eg + - CaC^ - c^Cc - CcCa is a ge- 
ometrical factor coming from the Cg orbital overlap, 
= —2t{ca ci) ^ Cc) ^ iJi is the bare dispersion in 
hole representation, and Cqc = cos /c^ with key are pro- 
jections of k on a, 6, c axes in the 2D hexagonal lattice 
of Co ions [Fig. [H^c)]. As the T-exciton has no disper- 
sion {eg-Cg hopping is zero in 90°-case), the selfenergy 
is momentum independent. Strong renormalization of 
the quasiparticle band and appearance of the incoherent 
sidebands as seen in Fig. [2] are the characteristic features 
of polaron formation. Excitations relevant here are the 
spin-state fluctuations of Co^+ ions, and a fermionic hole 
dressed by these excitations can be termed as spin-state 
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FIG. 3: (a) The physical picture behind Eq. (|3)). t2g electron 
of a Co^+ ion moves to eg level of the NN Co|+ ion (process 
1) and then to the t2g level of the next Co^^ neighbor (pro- 
cess 2). This is depicted in (b) as a motion of the hole-pair, 
(c) Singlet-pair motion via an intermediate state composed of 
S — 1 Co^^ ion (T-exciton) and triplet state of the two holes 
on (iA:)-bond. The relative amplitude A resulting from spin 
algebra is indicated, (d) The triplet-pair hopping amplitude 
A is three times smaller because of the destructive interference 
of contributions involving singlet and triplet (i/c)-bond. 



polaron. Physically, it is different from a typical mag- 
netic polaron formed in Mott insulators with a magneti- 
cally active ground state, while NaCo02 is a nonmagnetic 
band insulator. 

At large x ~ 1 limit, dilute polarons are readily 
trapped by a random potential of Na-vacancies^^^^^. 
When the binding is strong, physics is local and a polaron 
takes a form of hexagon-shaped S = 1/2 object where a 
hole is oscillating to optimize both t and i channels. Our 
model provides a microscopic basis for spin-polarons in- 
troduced on experimental grounds^^^ii^ and discussed in 
detail in Refs. I29II32I . When the density of polarons is 
increased (as x decreases), they start to overlap forming 
narrow bands. Eventually, the polaron picture breaks 
down and a correlated Fermi-liquid emerges when x is 
further reduced. 



IV. EFFECTIVE INTERACTION BETWEEN tag 
HOLES 



In the Fermi liquid regime, the Eliashberg- formalism, 
where the phonon shake-up processes (triggered by T- 
exciton) can also be incorporated, would be the best 
strategy. However, there are delicate constraints to han- 
dle: a lattice site cannot be occupied by two holes or 
by a hole and T-exciton simultaneously. For the sake 
of simplicity, we derive an effective fermionic interaction 
in a second order perturbation theory in t by consid- 
ering the local virtual process depicted in Fig. [3l(a,b). 



This way, all the constraints in the intermediate states 
are treated explicitly. Such a perturbative treatment is 
valid as long as a polaron binding energy (an energy 
gain due to the i process) is small compared to a bare 
bandwidth W 9t in a triangular lattice). From a 
self-consistent Born approximation discussed above, we 
obtained Ei, ^ 2M - 0.3VI^ for i = Et = 3t used in 
this paper. For this set of parameters, we can therefore 
integrate out a virtual spin states perturbatively. 

As a result, we arrive at the following effective Hamil- 
tonian in two equivalent forms: 



(a) 



local non-local 



(3) 



^eff = ^ COs{(l)ij - (t)jk) 

{ijk) 

= V ^ cos{(j)ij - (j)jk) [rijUik - \sjSik\ . (4) 



{ijk) 



We introduced here a constant V = P/Et- Sites 
i k are the nearest neighbors of site j. The an- 
gles (j) G (27r/3, 47r/3, 0) are selected by the orienta- 
tion of the bonds (ij) and {jk) as already explained. 
No-double-occupancy constraint on / is implied when 
using this effective Hamiltonian. The Hamiltonian in 
Eq. ([3|) describes the motion of the spin-singlet Sij = 

(Mfji - fiifj})/V2 and spin-triplet fij = {fi^fj^, . . .} 
Co^+-Co^+ pairs in a background of S = Co^+ ions, 
and may lead to the pairing instability as shown below. 

Alternatively, Eq. ^ represents the same interac- 
tion in a form of density-density and spin-spin correla- 
tions, emphasizing its relevance also to the charge and 
spin orderings. Note that Uik and Sik with i ^ k are 
the charge and spin densities residing on bonds, i.e., 

^ik = I Ea flfka^ ^ik = I Ea ^flfka (^hlle Tlj = Tljj 
and Sj — ^33 ^"^^ 

the usual on-site operators), so that the 
interaction acts between the local {on-site) and non-local 
(bond) operators. In a momentum space, Eq. (|4j) can be 
written as 



eff 



2V 



Q 



(5) 



with the operators = | Efc,a ^h+q^fl+q^afk^a^ = 

lT.k,a^^k+Q,kfUq^^U^^. etc. The formfactor Fj^^ j, = 
cos(A:a + k'^) + cos(A:5 + k^) + cos(A:c + k'J - Cac'^ - Chc'^ - 
Chc'^—Ccc'^—Ccc'^—Cac'^^ whcrc = COS /c^, originates from 
a peculiar bond-dependence of interactions in Eq. (|4]). It 
manifests again that the hq and Sq operators correspond 
to the particle-hole excitations that modulate the charge 
and spin bonds, respectively. 

To illustrate this unusual, nonlocal nature of correla- 
tions we show in Fig. [H the effect of the interaction on 
the spin susceptibility within the RPA approximation. 
The bare spin susceptibility [Fig. St^b)] is concentrated 
around the F point. When the interaction is switched on 
[Fig.HJc)], the 2kF ring in the susceptibility is enhanced. 
This suggests the fermionic 2 /cj^-inst abilities in a Fermi- 
liquid phase, consistent with a picture inferred from the 
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FIG. 4: (a) Diagrammatic representation of RPA equations 
for the spin susceptibilities involving the interaction Eq. JS} 
between Sq (local) and Sq (non-local) spin densities that re- 
side on sites and bonds, respectively. Bare and RPA-enhanced 
susceptibilities are represented by empty and shaded bub- 
bles respectively, (b) Map of bare x's 1^ for Ud — 0.5 (the 
average Co- valency 3.5) at T = 0.025t and uo = 0.005t. 
(c) Corresponding RPA-enhanced susceptibility calculated at 
P / Ft = 3t. The interaction enhances the susceptibility at 
the 2k F ring which (at given density na = 0.5) nearly matches 
the Brillouin zone boundary. 



experiments^. Interestingly, the RPA-spin susceptibility 
at Ud = 0.5 is most enhanced near the M point, i.e. 
near the observed magnetic Bragg peak position^, rather 
than at K typical for the AF Heisenberg spin system. 
In order to study the spin ordering at Ud = 0.5 more 
quantitatively, one should take into account also the Na 
ordering^ which breaks a hexagonal symmetry of the un- 
derlying Fermi-surface. 



V. SUPERCONDUCTIVITY DUE TO THE 
PAIR-HOPPING INTERACTION 

Now, we consider the above Hamiltonian in the context 
of superconductivity. It is evident from Eq. (|3]) that spin- 
singlet pairs gain much more kinetic energy than triplets. 
As explained in Fig. [3l this nontrivial result originates 
from a quantum interference between different realiza- 
tions of the virtual i process. The S = 1 T-excitation 
in the intermediate state is fully transparent for singlets 
which equally use all three Sz = ±1,0 states of T. How- 
ever, in the case of triplet pairs, there exist two quantum 
paths contributing with opposite signs, which results in a 
"spin blockade" for the motion of triplets. Alternatively, 
it can be said that the = Co^^ states move more 
coherently when the S = 1/2 background is in a singlet 
state. The difference from cuprates is that singlets are 
formed here not due to the superexchange (in cobalt ates, 
J is small^^) but because of the gain in the kinetic energy 
associated with i hoppings. 
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A mean-field BCS analysis of Eq. (|3]) shows that iJgff 
supports either extended s-wave singlet SC with the 
gap function (x j{k) = ^^/2/3{ca + C5 + Cc), or doubly- 
degenerate spin-triplet p-wave pairing with jx,y{k) = 
{{sa-Sb), {2sc-Sa-Sb)/V3}, where = sin/cc,. The d- 
wave channel is repulsive, while /-wave one is attractive 
but too weak in the physically reasonable doping range. 
We estimated the from 
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where Va is either VotV/3, and the corresponding form- 
factors are ^{k) or ^x,y{k) for the singlet 5- wave and 
triplet p-wave pairing, respectively. To account for the 
no-double-occupancy constraint, the fermionic dispersion 
as well as the pair-hopping amplitude are renormalized 
by the Gutzwiller factor^^ gt = 2nd/ (1 + ^d) as {i^V) = 
{gt(,^gtV)^ where rid is the relative fraction of Co^+ ions. 
(The reported Co-valences -3.4l^ -3.3^^, -3.46^^ op- 
timal for SC translate then to Ud = 0.6,0.7,0.54). In 
the momentum summation, we have introduced a cutoff 
equal to the excitation energy Et- 

We solved Eq. (|6]) at V = 3t (as it follows from 
i = Et = St used in previous sections). In terms 
of the BCS-coupling constant, this translates into A = 
VN = VN ~ 1/3 considering the density of states 
N ~ 1/W ~ l/9t. Therefore, the present formulation 
in terms of an effective fermionic Hamiltonian (3) should 
give a reasonable results. At larger values of we en- 
counter a strong coupling regime where one should use 
instead the original model (1) and treat a virtual spin 
states explicitely. This limit remains a challenging prob- 
lem for future study. 

The resulting Tc values from Eq. (|6|) are presented in 
Fig. [5](a) as solid lines. As expected, the highest Tc val- 
ues are found in the singlet channel, increasing with Co^+ 
density due to the formfactor effect, until SC disappears 
8it Ud = 1 limit. A weak triplet pairing is present thanks 
to its formfactor matching well the Fermi surface, but it 
is expected to be destroyed by {e.g. Na) disorder. (We 
should notice that these trends are based on the present 
mean- field decoupling which ignores a collective spin fluc- 
tuations. One can speculate, for instance, that the triplet 
pairing may be supported by a ferromagnetic fluctuations 
within the C0O2 planes observed^ at large Ud limit). 

As the SC pairing considered here is due to the pair- 
hopping. Coulomb repulsion between the holes will op- 
pose it. This is not a big trouble at high density of 
Co^+ spins (as they cannot avoid themselves) but be- 
comes a severe issue in a spin-diluted regime at large Ud, 
where Coulomb repulsion reduces the process described 
in Fig. [3] hence the amplitude V. Instead, the forma- 
tion of spatially separated spin-polarons (Fig. [2|) is fa- 
vored, and competing orderings take over, such as an in- 
plane ferromagnetism induced by a residual interactions 
between spin-polarons^^. To include the effects related 
to the Coulomb repulsion in the Gutzwiller fashion, we 
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FIG. 5: (a) Tc in the extended s-wave and p-wave chan- 
nels. The complete profile of the dominant, s-wave Tc curve 
is shown in the left inset together with 7^ (in arbitrary units) 
on the Fermi surface. The dashed {f3Vc — 1.5) and dotted 
{l3Vc — 3) Tc curves are calculated including NN Coulomb 
repulsion which reduces the pairing interaction V at large rid • 
Shaded regions indicate the observed competing orderings (in- 
cluding the spin-charge order at na = 0.5). (b) Probability 
ratio p{nd) (see text for definition) renormalizing the pair- 
ing interaction at different values of NN Coulomb repulsion 
relative to the effective temperature 1/(3 oc bandwidth. The 
feature at = 1/3 for large Vc manifests a honeycomb- 
lattice formation where each Co^^ (o) has the maximum pos- 
sible number of neighboring Co^^-Co^^ pairs (•-•). Above 
rid = 2/3, Co^^ holes can avoid each-other completely if Vc 
is sufficiently large. 



use an additional multiplicative factor reflecting the sup- 
pression of the probability Pijk of having the required 
Co^^ — Co^^— Co^^ configuration. We have determined 
this probability using a classical Monte-Carlo simulation 
of hardcore particles with NN Coulomb repulsion Vc- 
The simulations were performed at different "effective 
temperatures" 1//3 imitating the kinetic energy (of the 
order of bandwidth) which competes with the Coulomb 
repulsion in the real system. Plotted in Fig. [5fb) is 
the probability ratio p{nd) = PijkiVc) / PijkiVc = 0) 
for several values of pVc- The corresponding Tc curves 
calculated with V p{nd)V locate the SC-dome near 
the valence 3.4, in a remarkable correspondence with 
exp er iment— iiiii^ . 

Finally, our t — i model provides a clear hint on the 
role of water-intercalation needed for SC in Naa;Co02. 
Without water, a random Na-potential induces some 
amount of spin-polarons locally (the origin of "Curie- 
Weiss metal"-) which suppress the pairing among the 
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FIG. 6: (a) Band structure for different values of t' /t 
(rid = 0.6) and the corresponding effect on the s-wave transi- 
tion temperature. As the band dispersion near V point comes 
closer to the Fermi level, it can exploit the larger formfac- 
tor and Tc increases as shown in the inset. The increase is 
also partially induced by the decreased Fermi velocity along 
F-M direction, (b) Formfactor of the extended s-wave pairing 
interaction. The inset shows the contours for 7I = 1, 2, 3, 4, 5. 



remaining fermions the usual way. Once this potential 
is screened-out by the water layers, an intrinsic ground 
state of C0O2 planes as in Fig. [5] is revealed. (This in- 
terpretation of the water effect is consistent with the 
absence of superconductivity in the monolayer hydrate 
of Naa;Co02, where the water resides in the Na layers.) 
The remaining "enemy" of SC is the Coulomb repulsion 
which prevents the pairing of dilute Co^+ fermions and 
supports the formation of spin-polarons and magnetism 
instead. More pronounced polaron physics (because of 
the presence of large S = 1 T-exciton and narrow band- 
width) explains why Tc in cobaltates is low compared to 
cuprates. 

Another mechanism for the water effect is provided 
by the band-structure calculations^^ that indicate a sub- 
stantial flattening of the aig band-top and a reduction 
of the band splitting when the water-layers are present. 



To study the former effect, we include negative t' in our 
calculation. Due to the combined effect of better form- 
factor utilization in 5-wave channel and Fermi velocity 
reduction this enhances singlet pairing as presented in 
Fig. [6l In triplet channel, on the other hand, the t'-effect 
is weaker as it includes only the latter factor, i.e., the 
enhancement due to the reduced Fermi velocity (p-wave 
formfactor utilization is not very sensitive to t'). 



VI. CONCLUSIONS 

To summarize, we have presented t — i model for 
NaajCo02 which is based on the spin-state quasidegen- 
eracy of CoOe octahedral complex in oxides and the spe- 
cific lattice geometry of the C0O2 planes in layered cobal- 
tates. The model naturally explains the strong correla- 
tions found in the sodium-rich region due to the spin- 
polaron formation. We derived effective interactions in a 
Fermi-liquid regime and discussed their impact on spin 
fluctuations. The model predicts superconductivity me- 
diated by the spin-state fluctuations of Co^+ ions, at 
experimentally observed compositions. The basic idea 
behind the model is that due to the 90° d-p-d pathway 
in the edge-shared structure, the electron transport in 
NaajCo02 is entangled with low- lying S = 1 magnetic 
states of Co^+ ions which become an essential part of 
the Naa^Co02 physics. Given the simplicity and experi- 
mentally motivated design of the model, its success can 
hardly be accidental. Therefore, t — i Hamiltonian can be 
regarded as a basic minimal model for Naa;Co02. It may 
also have broader applications, e.g., in oxides of Rh and 
Ir ions with a similar spin-orbital structure and lattice 
geometry. 
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